####### EXERCCIOS 7 - ANA CAROLINA PAVAN ###########

#### Altura na Infncia e Na Vida Adulta ###
altura_jovem<-c(39,30,32,34,35,36,36,30)
altura_adulto<-c(71,63,63,67,68,68,70,64)
relao<-lm(altura_adulto~altura_jovem)
plot(altura_adulto~altura_jovem)
summary(relao)
# H uma relao significativa, corroborada pelo p<0.05
abline(relao)
esperados<-2*altura_jovem
esperados
abline(a=0, b=2, col=3)
confint(relao)
# Meus resultados no corroboram a hipotese do senso comum, pois o valor sugerido no est no intervalo de confiana.

##### Seriemas e Carcars #####
#1
AVES<-read.table("aves_cerrado.csv", header=TRUE,sep=";")
AVES
is.na(AVES)
AVES$urubu[is.na(AVES$urubu)] <- 0
AVES$carcara[is.na(AVES$carcara)] <- 0
AVES$seriema[is.na(AVES$seriema)] <- 0
AVES$fisionomia[AVES$fisionomia=="ce"] <- "Ce"
table(AVES)
str(AVES)
AVES$fisionomia <- factor(AVES$fisionomia, levels=c("CL","CC","Ce"))
str(AVES)
aves_ce<-AVES[AVES$fisionomia=="Ce",]
aves_ce
aves_cc<-AVES[AVES$fisionomia=="CC",]
aves_cc
aves_cl<-AVES[AVES$fisionomia=="CL",]
aves_cl
relao_ce<-lm(aves_ce$seriema~aves_ce$carcara)
summary(relao_ce)
relao_cc<-lm(aves_cc$seriema~aves_cc$carcara)
summary(relao_cc)
relao_cl<-lm(aves_cl$seriema~aves_cl$carcara)
summary(relao_cl)
### No h relao para nenhuma das fisionomias!
par(mfrow=c(2,2))
plot(aves_ce$seriema~aves_ce$carcara,ylab="seriemas",xlab="carcars", main="Ce")
abline(relao_ce)
plot(aves_cc$seriema~aves_cc$carcara,ylab="seriemas",xlab="carcars", main="Cc")
abline(relao_cc)
plot(aves_cl$seriema~aves_cl$carcara,ylab="seriemas",xlab="carcars", main="Cl")
abline(relao_cl)

#2
diferena = lm(seriema~carcara, data = AVES)
anova(diferena)
# No entendi direito como fazer, mas atraves desse resultado no h diferena

#### Resduos de Iris #####
#1
require(datasets)
head(iris)
iris_setosa<- iris[iris$Species=="setosa",]
str(iris_setosa)
setosa_sepala<-lm(iris_setosa$Sepal.Width~iris_setosa$Sepal.Length)
plot(iris_setosa$Sepal.Width~iris_setosa$Sepal.Length)
abline(setosa_sepala)
summary(setosa_sepala)
#2
setosa_petala<-lm(iris_setosa$Petal.Width~iris_setosa$Petal.Length)
plot(iris_setosa$Petal.Width~iris_setosa$Petal.Length)
abline(setosa_petala)
residuo<- residuals(setosa_petala)
residuo
setosa_petala2<-lm(iris_setosa$Petal.Width~residuo)
plot(iris_setosa$Petal.Width~residuo)
abline(setosa_petala2)
summary(setosa_petala)
summary(setosa_petala2)

### FIM ###






